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The fluctuation theorem (FT), the first derived conse- 
quence of the Chaotic Hypothesis (CH) of jjj], can be con- 
sidered as an extension to arbitrary forcing fields of the fluc- 
tuation dissipation theorem (FD) and the corresponding On- 
sager reciprocity (OR), in a class of reversible nonequilibrium 
statistical mechanical systems. 
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A typical system studied here will be N point particles 
subject to (a) mutual and external conservative forces 
with potential V(q\, . . . , (fjy), (b) external (non conser- 
vative) forces, forcing agents, {Fj}, j = 1, . . . , N, whose 
strength is measured by parameters {Gj}, j = 1, . . . , s, 
and (c) also to forces {fij}, j = 1, . . . N, generating con- 
straints that provides a model for the thermostatting 
mechanism that keeps the energy of the system from 
growing indefinitely (because of the continuing action of 
the forcing agents). 

An observable 0({q, q}) evolves in time under the time 
evolution St solving the equations of motion: 



ml = -d ql V(q/) + Fj({G}) + ^ 



(1) 



with m = particles mass, and <fj "thermostatting" forces 
assuring that the system reaches a (non equilibrium) sta- 
tionary state. The time evolution of the observable O on 
the motion with initial data x = (q,q) is the function 
t — > O(Stx) so that the motion statistics is the probabil- 
ity distribution fi + on the phase space C such that: 



lim — 



Q{S t x) dt 



0(y)n + (dy)= f (O). 



(2) 



for all data x S C except a set of zero measure with 
respect to the volume ~p on C. The distribution /i + is 
assumed to exist: a property called zero-th law, ^,|). 

The thermostatting mechanism will be described by 
force laws tpj which enforce the constraint that the ki- 
netic energy (or the total energy) of the particles, or 
of subgroups of the particles, remains constant, |Q. It 
is convenient also to imagine that the constraints keep 
the total kinetic energy bounded (hence phase space is 
bounded). The constraint forces {0j} will be supposed 
such that the system is reversible: this means that there 
will be a map i, defined on phase space, anticommuting 
with time evolution: i.e. Sfi = i S-t- 

Reversibility is a key assumption, p]J 



In M the divergence of the r.h.s. of Eq. (|l]) is a quan- 
tity — o~(x) defined on phase space that has been identified 
with the entropy production rate. 

The chaotic hypothesis, (CH), of p] implies a fluctu- 
ation theorem or (FT) which, 0), is a property of the 
fluctuations of the entropy production rate. Namely if 
we denote (<t)+ = J c a(y)n+(dy) the time average, over 
an infinite time interval by Eq. (0) , then the dimension- 
less finite time average p = p{x): 



T/2 



r/2 



a(S t x)dt d = (cr) + p 



(3) 



has a statistical distribution n T (p) with respect to the 
stationary state distribution /i + such that: 



t{ct) +P 



log 



Xt(-P) 
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(4) 



provided (of course) (<r}+ > 0. Following jlj a reversible 
system for which (<t)+ > will be called dissipative. 

Ruelle's 77-theorem, ||, states that (er) + > unless 
the stationary distribution /i + has the form p(x)jl (dx) . 

Hence we shall suppose that the system is dissipative 
when the forcing G does not vanish and, without real 
loss of generality, that o~(ix) = —o~(x) and that cr(x) = 
when the external forcing vanishes, writing: 



a( 2 ;)=^G I J«(x) + 0(G 2 ) 



(5) 



Assuming that the cr—cr correlations have a fast decay 
(a fact to be expected if the (CH) is accepted, ||) then by 
a result of Sinai |6) the entropy verifies a limit theorem; 



lim — log 7r r (p) 

T — >OG T 



'UP) 



(6) 



where £(p) is analytic for p in the interval [— p* , p*] within 
which it can vary (model dependent) The function 

£(p) can be conveniently computed because its transform 
A(/3) = linv^oo ilog/e^k-^+Tiv^dp can be ex- 
pressed by a cumulant expansion. Once \{0) is "known" 
then Q(p) is recovered via a Legendre transform; £(p) = 
ma X/3 (/3<a) + (p-l)-A(/3)), gg. 

By using the cumulant expansion for A(/3) we find that 
A(/3) = i/3 2 C 2 + ^t/3 3 C 3 + . . . where the coefficients C 3 
are (a(S tl >(5, 2 ■) . . . >(•))? dt, . . . if (. . .)£ 

denote the cumulants of the variables a(x). 
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In our case the cumulants of order j have size 0(G : '), 
by Eq. (§), so that: 



C(p) 



2C 2 



(p-l) 2 + 0((p-l) 3 (G 3 ) 



(7) 



(note the first term in r.h.s. giving the central limit theo- 
rem). Eq. (§), together with the (FT) (g), yields at fixed 
p our main relations: 



(*)+ = -C 2 + 0(G 3 ) 



(8) 



We define, [Q, the current Ji(x) = da i 'y{x) and the 
transport coefficients Lij — dj{Ji(x))+\G=o and we study 
Lij. 

To derive (FD), we first look at the r.h.s. of the first 
relation in Eq. (||) discarding 0(G 3 ): the r.h.s becomes 
a quadratic form in G with coefficients: 



dt({J?{Sr)j9{-)) + -{J?) + {J$) + ) 



(9) 



On the other hand the expansion of (cr}+ in the l.h.s. 
of Eq. rt§) to second order in G gives: 



(10) 



because the first order term vanishes (by Eq. (|^), or 
(§)). The r.h.s. of @ is the sum of |GjGj times 
djdj J a(x)/i+(dx) which equals the sum of the follow- 
ing three terms: the first is J didja(x)/j,+ (dx), the sec- 
ond is J diu{x)djpL + {dx) + (i «-> j) and the third is 
/ a(x)didjfi + (dx)), all evaluated at G = 0. The first 
addend is (by time reversal), the third addend is also 
(as a = at G = 0). Hence: 

d 4 d»+| G=0 = (di(j9) + + d l (j j >) + )\ G=0 (ii) 

and it is easy to check, again by using time reversal, that: 
flj(J?>+|G=o - dj(Ji)+\a=o = L iS (12) 

Therefore equating the r.h.s and the l.h.s. Eq. (||), 
as expressed respectively by Eq. (0) and (O) we get a 



formula for the matrix L ^f^^ giving (FD) (i.e. Green- 
Kubo's formula) at least if i = j. 

We want to show that the above ideas also suffice to 
prove (OR), i.e. Lij = Lji. The main remark is that 
we can extend the (FT) theorem to give properties of 
the joint distribution of the average of cr, (|^), and of 
the corresponding average of Gjdja. In fact define the 
dimensionless j-current q — q(x) as: 



,t/2 



def. 



-T/2 



Gjd J a(S t x)dt = J G j {d J a)+ q 



(13) 



where the factor Gj is there only to keep a and Gjdja 
with the same dimensions. 

Then if n T (p, q) is the joint probability of p, q the same 
proof of the (FT) in jjj yields also: 

lim^log^M) =1 (14) 
r^oo r{a) + p n T {-p,-q) 

and the limit theorem in (fy is extended, ||, to: 

C(P>«) = lim -log7r r (p,g) (15) 

r — >oo T 

We can compute £(p, q) in the same way as £(p) by 
considering first the transform A(/?i,/?2): 



l im I [ e TVi(r-i)W)++fh («-V^°)+)n T (r>,q)dpdq 

r—*oo t J 

(16) 

and then the Legendre trasform: 



t(p, q) = max (ft (p - l)(a)+ + ft (g - l)(G^-<7> + - 

P1.P2 



A(ft,ft)) 



(17) 



The function A(/3), (3 = (ft, ft), is evaluated by the 
cumulant expansion, as above, and one finds: 



\0) = \0,C0) + O(G 3 ) 



(18) 



where C is the 2x2 matrix of the second order cumulants. 
The coefficient C\\ is given by G2 appearing in (0); C22 is 
given by the same expression with a replaced by Gjdja 
while C12 is the mixed cumulant: 



((a(Sr) G J d, J a{-)) + - (a(S t -)) + {Gjdja{-))+)dt 

(19) 



Hence if w = ( , ^ll'^, I wcuei: 
{q-l)(G j d J a) + f 



(( P ,q) = l(C- 1 w,w)+0(G 3 ) 



(20) 



completely analogous to ([7]). But the (FT) in (|l4|), im- 
plies that £(p, q)~C(~~Pj ~Q) ls Q independent: this means, 
as it is immediate to check: 

(C- 1 ) 22 (G 3 d 3 a) + - (C- 1 ) 2 i(a)+ = + O(G 3 ) (21) 

which because of (||), and of (C _1 ) 2 2 = Cn/detC, be- 
comes the analogue of (§): 



(G J -a j -(7)+ = -C12 + 0(G 3 ) 



(22) 
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and, proceeding as in the derivation of (|9|) through (JE 
(i.e. expanding both sides of ( |2^ ) to first order in the 
Gi's and using (|T^)) we get that di(djO-) + is given by the 
integral in ([)[). This means that = Lji and the (FD) 
theorem follows toghetr with the (OR). 

Thus Eq. (§),(|22|) and the ensuing (FD), and (OR), are 
a consequence of (FT), ([I]), and of its (obvious) exten- 
sion, (pT[), in the limit G — > 0, when combined with the 
expansion (^) for entropy fluctuations. Those theorems 
and the fast decay of the o~o~ correlations are all natural 
consequences of (CH) for reversible statistical mechani- 
cal systems, which is the starting point of our consider- 
ations. Note that reversibility is here assumed both in 
equilibrium and in non equilibrium: this is a feature of 
gaussian thermostat models, 0, but by no means of all 
models, |,§. 

Of course while the (OR) and (FD) only hold around 
equilibrium, i.e. they are properties of G-derivatives 
evaluated at G = 0, and the expansion for X(f3) is a 
general consequence of the correlation decay, the (FT) 
also holds far from equilibrium, i.e. for large G and can 
be considered a generalization of the (OR) and (FD). 
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